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The Pierre Auger (cosmic ray) Observatory provides a laboratory for studying fundamental physics 
at energies far beyond those available at colliders. The Observatory is sensitive not only to hadrons 
and photons, but can in principle detect ultrahigh energy neutrinos in the cosmic radiation. Inter- 
estingly, it may be possible to uncover new physics by analyzing characteristics of the neutrino flux 
at the Earth. By comparing the rate for quasi-horizontal, deeply penetrating air showers triggered 
by all types of neutrinos, with the rate for slightly upgoing showers generated by Earth-skimming 
tau neutrinos, we determine the ratio of events which would need to be detected in order to signal 
the existence of new non-perturbative interactions beyond the TeV-scale in which the final state 
energy is dominated by the hadronic component. We use detailed Monte Carlo simulations to cal- 
culate the effects of interactions in the Earth and in the atmosphere. We find that observation of 
1 Earth-skimming and 10 quasi-horizontal events would exclude the standard model at the 99% 
confidence level. If new non-perturbative physics exists, a decade or so would be required to find it 
in the most optimistic case of a neutrino flux at the Waxman-Bahcall level and a neutrino-nucleon 
cross-section an order of magnitude above the standard model prediction. 



I. INTRODUCTION 

Ultrahigh energy cosmic neutrinos (UHEC^) are ex- 
pected to be produced in association with the observed 
ultrahigh energy (charged) cosmic rays (UHECR) , either 
at the same sites responsible for UHECR acceleration, or 
via interaction of the UHECR during propagation, par- 
ticularly with the cosmic microwave background (CMB). 
These neutrinos are unique probes of new physics as their 
interactions are uncluttered by the strong and electro- 
magnetic forces and, upon arrival at the Earth, they may 
experience collisions with center-of-mass energies up to 
y/s < 250 TcV. However, rates for new physics processes 
arc difficult to test since the flux of cosmic neutrinos is 
virtually unknown. Interestingly, it is possible in prin- 
ciple to disentangle the unknown flux and new physics 
processes by using multiple observables [lHl]. 

The Pierre Auger Observatory provides a promising 
way to detect UHEC^ by looking for deeply-developing, 
large zenith angle (> 75°) or "quasi-horizontal" air show- 
ers 0. At these large angles, hadron-induced showers 
traverse the equivalent of several times the depth of the 
vertical atmosphere and consequently their electromag- 
netic component is extinguished before reaching the de- 
tector. Only very high energy muons survive past about 2 
equivalent vertical atmospheres. Therefore, the shape of 
a hadron-induced shower front is very flat and prompt in 
time. In contrast, a neutrino shower exhibits the roughly 
same morphology as a vertical shower. It is therefore pos- 
sible to distinguish neutrino induced events from back- 
ground hadronic showers. Moreover, because of full fla- 
vor mixing, tau neutrinos are expected to be as abundant 



as other species in the cosmic flux. z/ T 's can interact in 
the Earth's crust, producing r leptons which may decay 
above to the Auger detectors; such events will be referred 
to as "Earth-skimming" events @, Q ■ 

Possible deviations of the neutrino-nucleon cross- 
section due to new non-perturbative interactions 1 can be 
uncovered at the Auger Observatory by combining infor- 
mation from Earth-skimming and quasi-horizontal show- 
ers. In particular, if an anomalously large rate is found 
for deeply developing quasi-horizontal showers, it may be 
ascribed either to an enhancement of the incoming neu- 
trino flux, or an enhancement in the neutrino-nucleon 
cross-section (assuming non-neutrino final states domi- 
nate). However, these possibilities can be distinguished 
by comparing the rates of Earth-skimming and quasi- 
horizontal events. For instance, an enhanced flux will in- 
crease both quasi-horizontal and Earth-skimming event 
rates, whereas an enhanced interaction cross-section will 
also increase the former but suppress the latter, because 
the hadronic decay products cannot escape the Earth's 
crust. Essentially this approach constitutes a straightfor- 
ward counting experiment, as the detailed shower prop- 
erties are not employed to search for the hypothesized 
new physics. 



Throughout this paper we use this term to describe neutrino in- 
teractions in which the final state energy is dominated by the 
hadronic component. We are not considering here new "pertur- 
bative" physics e.g. (softly broken) supcrsymmetry at the TeV 
scale which would have quite different signatures in cosmic ray 
showers. 
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In this paper, we compute how many Earth-skimming 
vs. quasi-horizontal showers one would have to collect at 
the Auger Observatory to convincingly demonstrate new 
non-perturbative physics in which the final state energy 
is dominated by the hadronic component. We show that 
even a small number of events could be sufficient to rule 
out the standard model (SM). Thus the expected low 
neutrino "luminosity" is not at all a show-stopper, and 
the Observatory has the potential to uncover new physics 
at scales exceeding those accessible to the LHC. 

In order to demonstrate this, we first compute accep- 
tances for Earth-skimming and quasi-horizontal events 
using detailed models of the terrain in the vicinity of 
the Observatory as well as detailed simulations of the re- 
sponse of the Surface Array to highly inclined air showers. 
We then perform a likelihood analysis to determine the 
event counts needed to exclude the SM at various con- 
fidence levels. The analysis includes systematic effects 
both from theoretical uncertainties in the (perturbativc 
QCD) SM cross-section Q and from uncertainties asso- 
ciated with the detector response. 

The outline of the paper is as follows. In Sec. [TT| we 
discuss some possible new physics scenarios which could 
manifest themselves as non-perturbative interactions at 
LHC energies and beyond. Next, in Sec. [TIT] we describe 
the detailed Monte Carlo studies of the acceptance for 
Earth-skimming and quasi-horizontal showers under the 
assumption of SM interactions, including systematic un- 
certainties. Finally in Sec. IIVI we perform the statistical 
analysis to ascertain the discovery reach of the Observa- 
tory. Our conclusions are collected in Sec. [V] 



II. NEW NON-PERTURBATIVE PHYSICS AT 
THE LHC ENERGY SCALE AND BEYOND 

The analysis techniques described herein constitute 
an entirely general approach to searching for non- 
perturbative interactions in which the final state is 
dominated by hadrons, without any dependence on 
what hypothetical mechanism might actually cause the 
'hadrophilia'. In order to illustrate some possible new 
physics signals which may be accessible using these tech- 
niques, we consider below two interesting possibilities. 



A. TeV-scale mass black holes 

In D-dimensional scenarios with large-compact-extra- 
dimensions (of common linear size 2irr c ) the Planck scale 
is related to the fundamental scale of gravity (Md) ac- 
cording to [T3| 



M, 



pi 



(1) 



— 1/2 

If Md > Mw = G F ~ 300 GcV, microscopic black 
holes (BH) can be produced gravitationally in particle 
collisions with center-of-mass energies s > 1 TeV fill ] . 



Subsequently a TeV-scale BH would promptly decay 
via thermal Hawking radiation (l2j into observable 
quanta [l]|. (For M D = 1 TeV, the lifetime of a BH 
of mass 10 TeV is less than 10 -25 s.) Since gravitational 
coupling is flavor blind, a BH emits all the ps 120 SM 
particle and anti-particle degrees of freedom with roughly 
equal probability. Accounting for color and spin, we ex- 
pect « 75% of the particles produced in BH evaporation 
to be quarks and gluons, « 10% charged leptons, « 5% 
photons or W/Z bosons, and rj 5% neutrinos. Thus, 
TeV BH production and evaporation constitutes a clear 
example of beyond-SM non-perturbative physics. 

Although such BH production cross-section ~ M^ 1 
is 5 orders of magnitude smaller than the QCD cross- 
section ~ Aq^ d , it was proposed [l4| that such BHs could 
be produced copiously at the LHC, and that these spec- 
tacular events could be easily filtered out of the QCD 
background. This is possible by triggering on BH events 
with prompt charged leptons and photons, each carrying 
hundreds of GeV of energy. 

Cosmic ray collisions, with center-of-mass energies 
ranging up to 10 5 GeV, certainly produce BHs if the 
LHC does. The question is, can they be detected? Most 
cosmic rays are protons or heavier nuclei, which collide 
with hadrons in the upper atmosphere, producing cas- 
cading showers which eventually reach the Earth's sur- 
face. At energies of interest, however, the cosmic ray 
luminosity (L ~ 7 X 10" 10 (£/PeV)~ 2 cm" 2 s _1 , tak- 
ing a single nuclcon in the atmosphere as a target and 
integrating over 27T sr), is about 50 orders of magnitude 
smaller than the LHC luminosity, thus making it futile to 
hunt for BHs in hadronic cosmic ray interactions. On the 
other hand, neutrino interaction lengths are far longer 
than the Earth's atmospheric depth, although they would 
be greatly reduced by the cross-section for BH produc- 
tion fl5| . Cosmic neutrinos therefore could produce BHs 
with roughly equal probability at any point in the atmo- 
sphere. As a result, the light descendants of the BH may 
initiate low-altitude, quasi-horizontal showers at rates 
significantly higher than SM predictions. 

Analytic and numerical studies have revealed that 
gravitational collapse tak es p lace at high energies and 
small impact parameters jl6l . Il7j . In the course of col- 
lapse, a certain amount of energy is radiated in gravita- 
tional waves, leaving a fraction y = Mbh/V^ available 
for Hawking evaporation. Here, Mbh is a lower bound 
on the final mass of the BH and y/J = 2xmNE v is the 
center-of-mass energy of the colliding particles, taken to 
be partons. This ratio depends on the impact parame- 
ter of the collision, as well as on the dimensionality of 
space-time. 

The inclusive production of BHs proceeds through dif- 
ferent final states for different classical impact parame- 
ters b (l7|. These final states are characterized by the 



fraction y(z) of the initial \f§ which is trapped within 
the horizon. Here, z = b/b max _, and 6 max = \fF r s is the 
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maximum impact parameter for collapse, where 

VI 2 B -V D - 7 >/ 2 r(^)" 



1 

Mb 



Mr 



D-2 



(2) 



is the radius of a D-dimcnsional Schwarzschild BH [18j |. 
and F is the form factor (l7j . 

The y dependance complicates the parton model calcu- 
lation, since the production of a BH of mass A/bh requires 
that s be Mg H /?/ 2 (z), thus requiring the lower cutoff on 
parton momentum fraction to be a function of impact 
parameter. Because of the complexity of the final state, 
we assume that amplitude intercfcrcncc effects can be 
ignored and we take the vN cross-section as an impact 
parameter-weighted average over partonic cross-sections, 
with the lower parton fractional momentum cutoff deter- 
mined by x m i n = Mg^f/Mo. This gives a lower bound 
X = (aWn-Mo) /[y (z)s] on the parton momentum frac- 
tion x. All in all, the vN — > BH cross-section reads [l9[ 

<t(vN^BR)= [ 2zdz [ dxFirr 2 V/*(a:,Q) , (3) 
Jo Jx i 

where i labels parton species and the /i(x, Q) are parton 
distribution functions (PDF). 

As an illustration, we consider the D = 10 string in- 
spired scenario. For, Md = 1 TeV, x m i n = 1, and 
primary neutrino energy E v = 10 10 GeV, we obtain 
a{uN ->• BH) ~ 2 x 10 6 pb pj|. This is about two orders 
of magnitude above SM predictions. The BH production 
cross-section by UHECV scales as 



a(uN -> BH) cx 



Ml 



D-2 
D-3 



(4) 



A further suppression arises if x min is increased. For pa- 
rameters in the semiclassical regime (x m i n ^ 3 [2(|) the 
BH cross-section becomes comparable to the SM cross- 
section at Mjj ~ 2 TeV; this determines the multidi- 
mensional Planck scale to which Auger may be sensi- 
tive. However, the LHC will also be sensitive to extra- 
dimensional effects at a similar scale. It is interesting to 
consider whether Auger may have access to new physics 
beyond the reach of the LHC and we now discuss such a 
possibility. 



B. Sphalerons 

In the electroweak theory, non-trivial fluctuations in 
SU (2) gauge fields generate an energy barrier interpo- 
lating between topologically distinct vacua [2l[. An in- 
dex theorem describing the fermion level crossings in the 
presence of these fluctuations reveals that neither baryon 
nor lepton number is conserved during the transition, but 
only the combination B — L. Inclusion of the Higgs field 
in the calculation modifies the original instanton con- 
figuration [22j . An important aspect of this modification 



(called the "sphaleron" ) is that it provides an explicit en- 
ergy scale of about 10 TeV for the height of the barrier. 
This barrier can be overcome through thermal transitions 
at high temperatures [23[ , providing an important input 
to any calculation of cosmological baryogenesis. More 
speculatively, it has been suggested 24[ that the topo- 
logical transition could take place in two particle colli- 
sions at very high energy. The anomalous electroweak 
contribution to the partonic process can be written as 



<7i(s) = 5.3 x 10 3 mb • e ~ (47r/Qw) Fw(e) 



(5) 



where aw — 1/30, the tunneling suppression exponent 
i*V(e) is sometimes called the "holy- grail function", and 
e = V~s/{^m w /a w ) ~ ^/30 TeV. Thus, it is even 
possible that at or above the sphaleron energy the cross- 
section could be of C(mb) [25|]. Of particular inter- 
est to cosmic ray physicists would be enhancement of 
the neutrino cross-section over the perturbative SM esti- 
mates, say by an order of magnitude in the energy range 
9.5 < \og w {E v /GeV) < 10.5. With the methods out- 
lined in this paper, this can be detected as an anomalous 
ratio of quasi-horizontal Earth-skimming showers to up- 
coming showers. 

It was shown [2|| that for the simple sphaleron 
configuration s-wave unitarity is violated for \/§ > 
4-7T Mw / &w ~ 36 TeV. For lower parton subencrgies, the 
cross-section is exponentially damped to values well be- 
low the perturbative SM electroweak value [25|, [2(| . On 
the other hand, at the higher parton subenergies, the 
cross-section may well be dominated by non-spherically 
symmetric classical field configurations (27| ■ If for V§ > 
36 TeV we saturate unitarity in each partial wave then 
this yields a geometric parton cross-section 7ri? 2 , where 
R is some average size of the classical configuration. As 
a fiducial value we take the core size of the Manton- 
Klinkhamer sphaleron, R ~ 4/Mw — 10 -15 cm. In this 
simplistic model, the vN cross-section is 

<L^(Eu)=*R 2 f £ f(x)dx, (6) 



11 partons 



where x E 



s min /s = {m) 2 /2ME u ~ 0.065. In the 



region 0.065 < x < 3(0.065) the PDF for the up and 
down quarks is well approximated by / ~ 0.5/x, so the 
expression for the cross-section becomes 



uN 
black disk 



{E v ) 



ttR 2 (0.5) (In 3) (2/2) 
1.5 x 10" 



- 30 cm 2 



(7) 



where the last factor of 2/2 takes into account the 
(mostly) 2 contributing quarks (it, d) in this range of x, 
and the condition that only the left-handed ones con- 
tribute to the scattering. This is about 80 times the SM 
cross-section [t| HH. Of course this calculation is very 
approximate and thecross-section can easily be smaller 
by a factor of 10 (e.g., if R is 1/3 of the fiducial value 
used) . 
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FIG. 1: Topography of the Auger site according to CGIAR- 
CSI data. The center of the map corresponds to the centre of 
the Auger array (latitude = 35.25° S, longitude = 69.25° W). 
The Auger position is marked by a circle. 



III. ACCEPTANCE AND SYSTEMATIC 
UNCERTAINTIES 

To calculate the acceptance we perform detailed Monte 
Carlo simulations. The incoming neutrinos are propa- 
gated through the Earth's crust, Andes mountains, and 
the atmosphere using an extended version [29| of the 
code ANIS do[. For fixed neutrino energies, 10 6 events 
are generated with zenith angles in the range 60° — 90° 
(down-going showers) and 90° — 95° (upgoing showers) 
and with azimuth angles in the range 0° — 360°. Ncuti- 
nos are propagated along their trajectories of length AL 
from the generation point on the top of the atmosphere 
to the detector in steps of AL/1000(> 6 km). At each 
step of propagation, the vN interaction probability 

P{E V ,E U 8) ~ N A o&(E v ) p(Z) AL , (8) 

is calculated using the cross-section (ogjjj (£?„)) estimates 
of Ref. @, where p(Z) is the local medium density, Ei the 
energy of the outgoing lepton, and Na — 6.022 x 10 23 g _1 . 
The outcoming particle spectrum from vN interactions 
is simulated with PYTHIA [3l[ and tau decays are sim- 
ulated using the package TAUOLA [HJ . 

The flux of otutgoing leptons as well as their energy 
and the decay vertex positions are calculated inside a 
defined detector volume. The geometrical size of the de- 
tector volume is set to 3000 x 10 km 3 and it includes the 
real shape of the Auger Observatory on the ground. A 
relief map of the Andes mountains was constructed ac- 
cording to a digital elevation data of the Consortium for 
Spatial Information (CGIAR-CSI) [H]. The map of the 
area around the Auger site is shown in Fig. [TJ 

The detection volume corresponds to the so called ac- 
tive volume in which potentially detectable neutrino in- 
teractions are simulated. For a given incoming neutrino 
with energy E v the active volume is defined by a par- 
ticular plane A gen and distance AL. The plane A gen is 



the cross-sectional area of the detector volume and it is 
used as a reference plane for the generation of incoming 
neutrinos. The area depends on the zenith angle 6 of the 
incoming neutrino. The distance AL is the multiple, n, 
of the average lepton range (L l (Ei)). 

Earth-skimming events occur in the Earth's crust, 
and so the relevant neutrinos and taus sample only the 
Earth's surface density, p s ps 2.65 g/cm 3 . At these ener- 
gies, the tau's propagation length is determined not by 
its decay length but by its energy loss. The energy loss 
per unit length of crossed matter is usually approximated 
by a linear equation (continuous energy loss approach). 
The t lepton loses energy in the Earth according to 

dE 

=-(a T + f3 T E T )p s , (9) 

where the factor a T parametrizes the ionisation losses and 
P T the energy losses through brcmsstrahlung, pair pro- 
duction, and hadronic interactions. For E v = 10 7 GeV, 
a T is negligible and (3 T ps 0.8 x 10~ 6 cm 2 /g (34|. Hadronic 
interactions (i.e., lepton- nucleus inelastic interactions 
dominated by small values of the squared momentum 
transfer Q 2 ) are responsible for the largest and the most 
uncertain contribution [35j. Such an uncertainty in j3 T 
dominates the systematic errors in the estimate of the 
neutrino event rates. 

To investigate the response of the Auger detector, we 
generate the lateral profiles of the shower development 
using the output of PYTHIA and/or TAUOLA as in- 
put for AIRES [3(|. The showers induced by the prod- 
ucts of up-going decaying tau leptons, with energies from 
0.1 EeV to 100 EeV and decay position at altitudes rang- 
ing from to 3500 m above sea leavel, are simulated 
in steps of 100 m. At each altitude 40 events are gen- 
erated to cover the tau decay channels implemented in 
ANIS 29|. In the case of down-going showers, the de- 
cay altitudes range from ground level up to the upper 
atmosphere. 

The response of the surface detector array is simu- 
lated in detail using the Offline simulation package [37l ]. 
Besides the standard procedure to simulate the spacial 
and temporal signal response we have added the simula- 
tion of atmospheric background muons to study the im- 
pact on the neutrino identification, since such accidental 
muons might be wrongly classified as shower particles. 
The background from hadronic showers above 10 8 GeV 
is estimated to be 0(1) in 20 years. At E v ~ 10 10 GeV 
the cosmic ray flux is ps 10 6 time smaller, so the ex- 
pected background for the energy bin considered here 
(9.5 < log 10 (£„/GeV) < 10.5) is negligible. 

The expected neutrino event rate (of flavor a) in the 
detector volume is found to be 

iVacc 

Nu a = F^^Pi, (10) 
where 7V acc is the number of events triggering the detector 
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and passing all quality cuts of the cascade analysis. Here, 

{• E m ax /""max 

F? = N gc i AT / (E u ) dE v / A gcn (9) dtt , 

(11) 

dfl is the solid angle, AT the observation time, N gen is 
the number of generated events from surface A gen , and 
we take the neutrino flux $>q™(E v ) to be isotropic. We 
further assume v e : : v T ~ 1 : 1 : 1, which is gener- 
ally thought to be the case if the neutrinos are produced 
predominantly thorugh pion decay. In order to ascertain 
the systematic uncertainties associated with our lack of 
knowlcgc of the dependence of the flux on energy, we con- 
sider three scenarios which plausibly bracket the range of 
possibilities: 

1. %«{E v ) = {C/E Q )E-\ 

2. ^(E U )=CE- 2 , 

3. $ V "(E V ) = (C/E )E- 3 , 

4. %"{E V ) = CE- 2 eK P [-log 10 (E 1/ /E f/(2a 2 )}, 

where C = 2.33 x 1CT 8 GeVs^ 1 cm~ 2 sr" 1 , E Q = 
10 10 GcV, a = 0.5 GcV. This normalization (2) con- 
stitutes the common benchmark, the so-called 'Waxman- 
Bahcall bound' [H . The factor F™ of Eq.[TT]is chosen to 
yield the total number of events per year. The expected 
rates for the entire range over which Auger is sensitive 
are given in Table U and the rates for the high energy bin 
considered in the following study are given in Table HU 



TABLE I: Expected events per year (Ni) at Auger in the 
energy range 8 < log 10 (_E„/GeV), for various incident zenith 
angle (9) ranges, assuming the Waxman-Bahcall flux. 



flux 


up-going 


down-going 


ratio 


(2) 
(2) 


e N Vt 


e N Ve N Vt N Uu N vm 




90-95 0.68 
90-95 0.68 


60-90 0.134 0.109 0.019 0.262 
75-90 0.075 0.071 0.011 0.157 


2.58 
4.27 



TABLE II: Expected events per year (Ni) at Auger in the en- 
ergy range 9.5 < log 10 (i?„/GeV) < 10.5, for various incident 
zenith angle (9) ranges and the 4 flux models considered. 



flux 


up-going 


down-going 


ratio 




9 N Vt 


e 


N Ve 


N Vt 


Nu„ 




Nr/N^ u 


(1) 


90-95 


0.14 


60-90 


0.059 


0.049 


0.011 


0.12 


1.14 


(2) 


90-95 


0.15 


60-90 


0.059 


0.049 


0.096 


0.11 


1.33 


(3) 


90-95 


0.23 


60-90 


0.079 


0.062 


0.0123 


0.15 


1.53 


(4) 


90-95 


0.12 


60-90 


0.046 


0.037 0.0080 


0.091 


1.33 


(1) 


90-95 


0.14 


75-90 


0.027 


0.031 


0.0056 


0.064 


2.14 


(2) 


90-95 


0.15 


75-90 


0.026 


0.029 


0.0048 


0.060 


2.47 


(3) 


90-95 


0.23 


75-90 


0.036 


0.041 


0.0062 


0.083 


2.75 


(4) 


90-95 


0.12 


75-90 


0.021 


0.024 


0.0040 


0.049 


2.45 



Hereafter we consider ^ a (E v ) cx E~ 2 as our nom- 
inal spectrum. We then estimate systematic uncer- 
tainties associated with: different assumptions of the 



spectrum sha pe, different parton distribution functions 
(GRV92NL0J3| and CTEQ66c g3), and different esti- 
mates on j3 T (35( . The contribution of different systematic 
errors are listed in Table IIIII 

TABLE III: Contributions to the systematic uncertainty on 
the Earth-skimming to quasi-horizontal event ratio. We have 
considered the energy range 9.5 < log 10 (_E„/GeV) < 10.5 and 
the zenith angle range 75° < 9 < 90°. 



ratio 


flux 


PDF 


fir 


sum 




+ 11% 


0% 


+24% 


+ 26% 


2.47 








2.47 




-13% 


-21% 


-25% 


-35% 



IV. AUGER DISCOVERY REACH 

Consider a flux of neutrinos with energy in the range 
10 9 - 5 GeV < E v < 10 10 - 5 GcV. In the SM, the interac- 
tion path length is 

L v cc = [Nap^cc)- 1 , (12) 

where (Tq C is the charged current cross-section for E v = 
E . (We neglect neutral current interactions, which at 
these energies serve only to reduce the neutrino energy 
by approximately 20%, which is within the systematic 
uncertainty.) For E a - 10 10 GeV, L v cc - 0(100) km. 
Supplemented by the possibility of new non-perturbative 
physics, the interaction path length is 

L^ ot = [N A p s (^ cc +<p)r 1 , (13) 

where <r^ P is the new physics contribution to the cross- 
section for E v = Eq. 

The maximal path length for a detectable r is given 

by 

L T = -t— — In (E max / E m [ n ) , (14) 

PrPs 

where E max sa Eq is the energy at which the tau is cre- 
ated, and E m { n is the minimal energy at which a r can 
be detected. For E ma ^/E m [ n = 10, L T = 11 km. 

Given an isotropic v r + v T flux, the number of taus 
that emerge from the Earth with sufficient energy to be 
detected is proportional to an "effective solid angle" 

f2 eff = J dcos9d(j) cos 9 P(9,(j)) , (15) 

where 

P(6, 4>) = / ^e" z/L - Q[z-(t- L T )\ (16) 

JO L GC 

is the probability for a neutrino with incident nadir an- 
gle 9 and azimuthal angle <j> to emerge as a detectable 
r. (In Eq. tfrB")) , for the reasons noted above, we have 
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FIG. 2: Projected determination of neutrino fluxes and cross-sections at yfs ss 250 TeV from future Auger data. The different 
shaded regions indicate the 90%, 95%, 99% and 3<r confidence level contours in the $" a11 /^q" 11 — ctnp/o-cc plane, for Ags = lj 
Aqh = 10 (left), iVes 8 = 1, ^qh = 7 (middle), and N£f = 1, = 5 (right). The dashed line indicates the result of 

including the systematic uncertainty on the NLO QCD CC neutrino- nucleon cross-section [42l |. 



neglected the possibility of detectable signals from new 
non-perturbative physics by Earth-skimming neutrinos.) 
Here £ = 2i? ffi cos 9 is the chord length of the intersec- 
tion of the neutrino's trajectory with the Earth, with 
R® ~ 6371 km the Earth's radius. Evaluating the inte- 
grals, we find [l| 



f2 e ff = 27T- 



T '■ 



J cc 



- 1 



T '■ 



2Ra 




-2R 9 /L» 



(17) 



At the relevant energies, the neutrino interaction length 
satisfies L" ot <C R®. In addition, for L" ot ^> L T . valid 
when the cross-section enhancement is significant but not 
so large as typical hadronic cross-section, Eq. (fTT| simpli- 
fies to @ 



fi 



eff 



2tt 



r i/2 tt 

4R^L£ C 



(18) 



Equation (|18|) gives the functional dependence of the 
Earth-skimming event rate on the new physics cross- 
section. This rate is, of course, also proportional to the 
source neutrino flux &" lx11 at Eq. Given these inputs, 



N E s ~ Ces -j-T^n 



'CC 



(19) 



where Ces = 0.15 is the number of Earth-skimming 
events expected for a fiducial flux $Q aI1 in the abscence 
of new physics. 

In contrast to Eq. ([T9"l) . the rate for quasi- horizontal 
showers has the form 



QH 



V <J CC (7 N p 

Or" 



(20) 



where Cqh = 0.06 for the Auger Surface Array, as deter- 
mined in Sec. IIIII 

Given a flux $ lya11 and new physics cross-section cr^ P , 
both N-es an d A^qh are determined. On the other hand, 
given just a quasi-horizontal event rate Aqh, it is im- 
possible to differentiate between an enhancement of the 
cross-section due to new physics and an increase on the 
flux. However, in the region where significant event rates 
are expected the contours of A^qh and A^es, given by 
Eqs. (fTTij) and ([2"U|) . are more or less orthogonal and pro- 
vide complementary information. With measurements of 
Ng£ and Agjf , both cr^ p and <J>^" may be determined 
independently, and neutrino interactions beyond the SM 
may be unambiguously identified. 

We now turn to determining the projected sensitiv- 
ity of Auger to neutrino fluxes and cross-sections. The 
quantities Aes and Aqh as defined in Eqs. (fT9)l and ([20]) 
can be regarded as the theoretical values of these events, 
corresponding to different points in the (j)"*" /$g aI1 — 
cnp/ccc parameter space. For a given set of observed 
rates A$f and N$$, two curves are obtained in the two- 
dimensional parameter space by setting N^ s = A^es and 
^Qjf = ^QH- These curves intersect at a point, yielding 
the most probable values of flux and cross section for the 
given observations. Fluctuations about this point define 
contours of constant \ 2 i n & n approximation to a multi- 
Poisson likelihood analysis. The contours are defined by 



X 



E 2 i N * 



N 



obsl 



+ 2AT° bs In [N° hs /Ni] , (21) 



'cc 



where i = ES, QH [41(. In Fig. [2j we show results 
for three representative cases. Assuming (A^o 8 = 1, 
N$£ = 10), (Nvf = 1, = 7), and (JVg§« = 1, 

N^g = 5) we show the 90%, 95%, 99% and 3a CL con- 
tours for 2 d.o.f. (x 2 = 4.61, 5.99, 9.21, and 11.83, re- 
spectively). For N^ s = 1 and = 10, the possibility 
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of a SM interpretation along the a^ P = axis (taking 
into account systematic uncertainties) would be excluded 
at greater than 99% CL for any assumed flux. The power 
of the Earth-skimming information is such that the best 
fit consistent with the SM would require a flux of about 
50 times the Waxman-Bahcall flux, which is already ex- 
cluded by present limits |43| . 



cross-section about an order of magnitude above the stan- 
dard model prediction). The proposed Northern Auger 
site [44[ (which has not been optimized for neutrino stud- 
ies) would reduce this time by about a factor of 2. Any 
hint of such an important signal would provide an impe- 
tus to infill the array to increase the neutrino acceptance. 



V. SUMMARY 

We have re-examined a technique to search for new 
physics at sub-fermi distances. The strategy involves 
determining the ratio of quasi-horizontal to Earth- 
skimming showers initiated by cosmic neutrinos which 
would need to be detected by the Pierre Auger Obser- 
vatory in order to signal the existence of exotic non- 
perturbative interactions beyond the TcV-scale. We per- 
form Monte Carlo simulations of neutrino interactions in 
the Earth and in the atmosphere, and realistic simula- 
tion of the detector acceptance using the Auger Off line 
software. We find that observation of 1 Earth-skimming 
and 10 quasi-horizontal events would exclude the stan- 
dard model at the 99% confidence level. If new non- 
perturbative physics exists, a decade or so would be re- 
quired to uncover it in the most optimistic case (cos- 
mic neutrino flux at the Waxman-Bahcall level and vN 
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